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ABSTRACT 

We employ a suite of 75 simulations of galaxies in idealised major mergers (stellar 
mass ratio ^2.5:1), with a wide range of orbital parameters, to investigate the spa¬ 
tial extent of interaction-induced star formation. Although the total star formation in 
galaxy encounters is generally elevated relative to isolated galaxies, we find that this 
elevation is a combination of intense enhancements within the central kpc and moder¬ 
ately suppressed activity at larger galacto-centric radii. The radial dependence of the 
star formation enhancement is stronger in the less massive galaxy than in the primary, 
and is also more pronounced in mergers of more closely aligned disc spin orientations. 
Conversely, these trends are almost entirely independent of the encounter’s impact pa¬ 
rameter and orbital eccentricity. Our predictions of the radial dependence of triggered 
star formation, and specifically the suppression of star formation beyond kpc-scales, 
will be testable with the next generation of integral-field spectroscopic surveys. 
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panions IIBarton et al. 2000l: Barton Gillespie et al.l 

20031; 

Galaxy encounters are recognized as a leading mechanism 

Lambas et al.l l2003l: lAlonso et al.l l2006l: IWoods et al. 

2006|; 

for triggering star formation (ISanders & Mirabel! Il996f). 

Ellison et al.ll2008l: Patton et al.ll201l]: IScudder et al 

2 OI 2 I; 

Numerical simulations support this picture: interaction- 

Scott & Kavirail 2014h. Relevant related works include in- 


driven non-axisymmetric gravitational torques funnel co¬ 
pious amounts of gas into the centr al regions, fuel¬ 
ing powerful bursts of s t ar fo rmation (iHernauistI 1 19891 : 
iBarnes fc HernauistI 1 19911 . Il996ll . However, this does not 
necessarily mean that starbursts are ubiquitous in interact¬ 
ing galaxies. Triggering depends on many factors, i nclud¬ 
ing the specihc mergi n g geo r netrv (ICox et al.l 1200^ . l2008l : 
IPi Matteo et ah] 120071 . koOSi : iTorrev et ajj _ 20U ) and the 


prop e rties of the progenitor galaxies ( iMihos fc HernauistI 


1 19961 : ISpringell l2000l : ISpringel fc HernauistI 2005ll . Even if 
triggering occurs, the observability time-scale of the burst 
may be shorter than the duration of the interaction. In other 
words, a galaxy in an encounter may be identified as normal 
(or even passive) if observed long after the peak of activity. 

Despite these issues, surveys detect (on average) ele- 


vestigations on the role of environment ( 

Alonso et al. 20041: 

Ellison et al.l I 2 OIOI: lAlonso et al.l 12012 

: Kamoczvk et al.l 

2 OI 3 I: Ellison et al.ll2013f) - and extensions to wider separa- 

tions lout to ~150 kpc:|Patton et al.ll2013f) and high redshift 

llLin et al.ll2007l: iFreedman Woods et al.l 

2 OIOI: iHwang et al.l 

201l]:IWong et al. 2011): Kamoczvk et al 

2013h. 


To shed light on how merger-induced star formation 
unfolds, it is interesting to ask where the transformation of 
gas into stars is most efficient. Several works address this 
question by focusing on single individual systems in detail 
(e.g., the Antennae - see Section EJ. Unfortunately, large 
surveys are not yet able to provide a statistical view of the 
spatial distribution of star formation in galax ies with close 
compa nions. A crude attempt is performed bv lEllison et al.l 
(l2013ll . who measure the star formation rate (SFR) in aper- 
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2 J. Moreno et al. 


tures (‘fibres’) centred on galaxies, an d compare it to the 
total SFR of those galaxies (see also IPatton et alJ I 2 OIII : 
IScudder et alJEoid : IPatton et al.ll201.‘ih . They show that fi¬ 
bre SFR is more elevated (relative to control galaxies with¬ 
out close companions) than total SFR - suggesting that 
centrally-concentrated star formation in interactions is com¬ 
mon. Unfortunately, using a fixed angular fibre means that 
different portions of the target galaxies are covered - i.e., 
coverage depends on how large galaxies appear on the sky 
(i.e., their distance to us). This obstacle severely limits our 
ability to properly quantify the spatial concentration of star 
formation in interacting galaxies. 

Ideally, surveys equipped with the ability to map the 
location and kinematics of star forming regio ns would pro¬ 
vide b e tter clues on how th is pr ocess unfold s. Kewlev et al.l 
ll2010tl . iRupke et al. Idioia) and iRosa et ^ (120141 ') use HII- 
regions to find that interacting galaxies have shallow metal- 
licity gradients, which is consistent wit h inflow of gas tha t ul¬ 
timately ignites starbursts. Similarly, I Rich et ^ ll2012li use 
integral field spectroscopy (IFS) to create detailed metal- 
licity maps in interacting syste ms. More direct approaches 
include fKnapen fc Jamea ll2009ll . wh o use Ha imaging to es - 
timate SFR profiles in mergers - and iBellocchi et al.l ll2013ll . 
who employ IFS to determine changes in the 2D kinematic 
behaviour of ionised gas (Ha) as the merging sequence ad¬ 
vances. Their main limitation, however, is that the samples 
considered often contain too few galaxies. In other words, it 
is not clear that the information we infer from just a handful 
of c ases is universa l. _ 

ISchmidt et al.l ll2013li are the first to attempt using a rel¬ 
atively large sample: 60 mergers at 2 ; ~ 1.5, visually-selected 
with 3d-hst. These authors use near-infrared slitless spec¬ 
troscopy to measure the spatial extent of star formation 
(via Ha and [OHI] emission-line maps). Unfortunately, their 
scheme is very crude: it only checks if star formation occurs 
in a single galaxy, in both, or in the region connecting the 
two galaxies. In this paper, we underscore the need for obser¬ 
vations capable of measuring resolved star-formation maps 
in interacting galaxies. 

In this direction, the emergence of large multi-Integral- 
Field-Unit (IFU) spectroscopic progr ammes holds consid - 
erable promise. Surveys lik e CALIF^J (lSanchez_jt_alJ 20lJ), 
SAm| 3 (ICroom et al.l f20l3) . Ma NGifI llBundv et al.l 2015fi . 
and the future HECTOR survey (lUawrence et al. 2012f) will 
soon be able to analyse large samples of interacting galaxies 
with exquisite spatial detail. Indeed, both CALIFA (Barrera- 
Ballesteros et al., in prep) and SAMI are already on their way 
to analysing their respective samples of interacting galaxies 
(Jorge Barreda-Ballesteros and Iraklis Konstantopoulos, pri¬ 
vate communications). Beyond spect roscopic mapping^ mid- 
infrared imaging surveys like s4c 3 (ISheth et al.ll201(J ~) will 
also provide clues on the me rging sequence, and its impact 
on the structure of galaxies dKnapen et al.ll2014^ . 


^ The Calar Alto Legacy Integral Field spectroscopy Area survey 
(http://califa.caha.es ) 

The Sydney-Australia-Astronomical-Observatory Muti-object 
Integral-Field Spectrograph (http://saini-surv6y.org) 

® Mapping Nearby Galaxies at the Apache Point Observatory 
(https://www.sdss3.org/future/manga.php ) 

Spitzer Survey of Stellar Structure in Galaxies 
(http://www.cv.nrao.edu/~ksheth/S4G) 


With these surveys in sight, the time is right to con¬ 
duct resolved spatial studies of star formation in interacting 
galaxies with numerical simulations. The aim of this paper 
is to answer the following question: is star formation in in¬ 
teracting galaxies nuclear or extended? Also, which orbital 
parameters govern the spatial extent of star formation? 

This paper is organised as follows. We present our meth¬ 
ods in Section [ 2 ] Section [3] describes a case study, and Sec¬ 
tion |3] generalises to other merger configurations. We discuss 
our findings in Section [S] and summarise in Section!^ 


2 METHODS & DEFINITIONS 
2.1 The Model 


We use the smoothed-par ticle hydrodynamics (SPH) code 
gadget-3 (ISpringell l2005l i to run idealised galaxy merger 
simulations. Th ese include gravi ty, hydrodynamics, radia¬ 
tive gas cooling ( Katz_et_^ ]ll996ll. star form ation with asso¬ 
ciated feedback ( Springel fc Hern auistll20 03li. and supermas - 
sive black hole growth and feedback ( Di Matteo et al.ll2005l i . 
This model employs an equation of state parameter oi q = 
0.3 to handle the pressurization of de nse, star forming gas. 
Rega rding our formulation of SPH (ISpringel fc HernauistI 
l2002li . we do not expect our results to be se nsitive to de¬ 
tails of the hydro solver (iHavward et al.ll2014f) . 

This model has the advantage of being well numerically 
converged. By suppressing gas fragmentation, it ensures 
tha t star formation p roceeds at an efficiency consistent with 
the iKennicuttI (1 19981 ') (KS) relation. Therefore, the results 
presented in this paper are mostly dependent on our ability 
to resolve the mechanisms that redistribute the gas through¬ 
out the galaxy (e.g., bars, arms, tidal forces, gas shock heat¬ 
ing, etc.) with the local SFRs then being determined by an 
enforced volumetric KS relation. _ 

Other models in the literature (e.g ., ^^ssie^jt^y|20^; 
iPowell et al.ll20l3 : iRenaud et al.l[2014 : iRenaud et al.ji2015f ) 
attempt to determine the SFR efficiency in mergers by re¬ 
solving parsec scale physical processes. Here, we only hope 
to understand where the gas has moved during the merger, 
and therefore where we would expect SFR to occur. It is 
therefore worth cautioning that if our adopted KS relation 
is invalid in merging systems, then this could impact our 
results. However, a significant violation of the KS relation 
would be required to change our qualitative conclusions. 


2.2 Galaxy Merger Simulations 


In this paper, w e empl oy the simulation suite first discussed 
in IPatton et al.l (|2013l). Th is mer ger suite differs from that 
presented in Torrev et al.l ll2012ll only in the adopted ini¬ 
tial conditions and orbital parameters. The same simulation 
code and physics modules are adopted in both. As such, we 
direct the reader to lTorrev et al.l ll2012ll for a more complete 
description of our employed physics modules. 

We focus exclusively on mergers with stellar mass ratios 
~2.5:1 - where, initially, the larger (primary) galaxy has 
stellar mass of M* = 1.4 x 10^°M© and the smaller (sec¬ 
ondary) galaxy has a stellar mass of M* = 5.7 x 10®M©. 
These specific choices represent typical values in the Sloan 
Digital Sky Survey galaxy-pair catalogue of IPatton et al.l 
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Orientation 

Identifier 

4^1 

[degrees] 

Oi 

[degrees] 

4>2 

[degrees] 

O 2 

[degrees] 

e 

60 

30 

45 

-30 

f 

60 

60 

0 

150 

k 

-30 

-109 

-30 

71 


Tabl e 1. The merger orienta tions considered in this paper, drawn 
from iRobertson et alj ll2006l~l , and described in Figure ^ 


ll2013l) - and are also commonplace i n the cosmological 
galaxy-pair catalogue of iMoren o et al.| (|2013t), drawn from 
the Millennium Simulation I Springel et al .I I 2 OO 51 J 1 . 

For both galaxies, we adopt an initial stellar bulge- 
to-dis c ratio of Mbuige/Mdisc = 0.24, as in IPatton et al.l 
(l2013tl . The ini tial galax ie s are set up following the an¬ 
alytic wo rk of Mo et aP (Il998l ~l. via the procedure out- 
li ned in Springel et al.l (l2005l l. Simulated bulges follow 
a lHernauisd ( lOToh profile - and simulated discs have scale 
lengths of 2.1 kpc (for the more massive galaxy) and 1.5 
kpc (for the less massive galaxy), for both the gaseous 
and stellar components. The initial gas fraction is set at 
■feas = Afgas/Mdisc = 0.25, consistent with observations 
(ICatinella et ak Iffl). 


Each run has ~2.5xl0° baryon particles, yielding a 
baryon mass resolution of Mb 10“*, with a Plummer equiv¬ 
alent gravitational softening length of 50 pc. Each simula¬ 
tion has 2x10^ dark matter particles, yielding a dark matter 
mass resolution of Mdm = 2.5 x 10® Mq, with a Plummer 
equivalent gravitational softening length of 200. 

Our merger simulations are initiated by placing two 
otherwise stable galaxies on an interacting orbit. We con¬ 
sider 75 merger simulations in total in this paper, con¬ 
sisting of 25 different orbital configurations (variations in 
the orbital energy and angular momentum) with three dif¬ 
ferent alignments of the galaxy’s angular momentum rel¬ 
ative to the plane of the merger (see below). The 25 or¬ 
bital configurations are built by considering orbital ec¬ 
centricities t = {0.85,0.90,0.95,1.0,1.05} and Keplerian- 
inferred impact parameters of 6 = {2,4,8,12,16} kpc. 
These choices are consisten t with cosmological simulations 
dKhochfar fc Burkertll2006l ~l. 

We employ the following three merger orientations: 


these are the “e”, “f ’, and “k” orientations drawn 
from IRobertson et al.l (l2006fl . and summarised in Table [T] 
These orientations are selected to represent two strongly 
aligned discs (“e”-orientation), two nearly perpendicular 
discs (“f’-orientation), and two nearly anti-aligned discs 
(“k”-orientatio n). See Figured] fo r a schematic description 
(adapted from ll^orrev et ^l2^2ll . 



Figure 1. Schematic representation of angles defining the relative 
orientation of our merging galaxies. See Table [T] fo r the specific 
angles employed in this paper. Figure adapted from iTorrev et aU 
ll2ni2ll - see their Figure 6. ©AAS. Reproduced with permission. 


procedure becomes particularly tricky when the separation 
between the two galaxies is smaller than their typical sizes. 

Our approach is slightly different. For each galaxy, we 
only focus on the SPH particles contained within a sphere 
of 10-kpc radius centred around the corresponding super- 
massive black hole. We adopt this radius because our two 
galaxies in isolation form their stars entirely within 10-kpc. 
We also checked that the fraction of star formation taking 
place outside such spheres in the interacting case is neg¬ 
ligible for all of our runs. This approach is powerful be¬ 
cause it allows us to compare galaxies with companions to 
their isolated equivalents directly, on a region-by-region ba¬ 
sis. Breakdown occurs when the two spheres overlap: SPH 
particles are assigned to two galaxies simultaneously, lead¬ 
ing to double counting. To avoid this complication, we ignore 
those few snapshots where the separation between the black 
holes is less than 20 kpc. 

In this paper, we are primarily interested in the stages 
of merging where the two distinct galaxies can be identi¬ 
fied unequivocally. For this reason, we focus exclusively on 
the interacting phase: the period between first and sec¬ 
ond pericentric passage (with the provision that separation 
is greater than 20 kpc). This adopted approach facilitates 
comparing the spatial extent of star formation in interact¬ 
ing galaxies to that in their isolated counterparts. 


2.3 Galaxy Membership & the Interacting Phase 

The central goal of this paper is to map star formation in in¬ 
teracting galaxies, and compare to equivalent maps of their 
isolated counterparts. A nuisance of this procedure is how 
to assign SPH particles to each of the two galaxies. Some 
works assign membership by seeking the nearest supermas- 
sive black hole, which is treated as a proxy for the ce ntre of 
its host galaxy llTorrev et al.l[20l3 : |Patton et al.ll2013li . This 


3 RESULTS: A CASE STUDY 

This section describes a case study in full detail (eccentric¬ 
ity e = 1.05, impact parameter 6 = 16 kpc). This choice is 
not meant to be average. Instead, our goal is to maximise 
the duration of the interacting phase (prior to coalescence). 
All other runs exhibit qualitatively similar features, except 
that they are always interrupted by merging at an earlier 
time. We first focus primarily on the “e” orientation. See 
Section 143] for other orientations. 
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Figure 2. Density maps: gas [top and middle, blue) and star-forming gas [bottom, pink). Upper row displays 120 kpc X 120 kpc stamps 
of the interacting system. Middle and lower rows are centred on the secondary (smaller) galaxy, displayed on 20 kpc X 20 kpc stamps. 
Columns (left-to-right) show the following stages: incoming (a); first passage (b); apocentre (c); and second approach (d). First passage 
produces tidal tails, followed by a nuclear starburst and a mild star-forming contribution at larger galacto-centric radii. 


3.1 Mapping Star Formation 

Figure [2] shows density maps of the gas (blue) and star form¬ 
ing gas (pink) for our case study. The upper row shows 120 
kpc X 120 kpc stamps, depicting the evolution of the in¬ 
teraction on extra galactic scales. The middle and bottom 
rows show 20 kpc x 20 kpc stamps centred on the secondary 
(smaller) galaxy. Star-forming gas (bottom row) only traces 
the densest gas (middle row, in blue), as expected in ou r 
iKennicut^ ll 19981) based model llSpringel fc Hernauisdl200^ . 

Columns (a)-(d) represent various stages of interaction: 

• (a) Incoming Phase: Before the interaction. Left 
alone, these galaxies retain their original morphology, and 
exhibit declining star formation as a function of time. 


• (b) First Pericentric Passage: The galaxies exhibit 
short-lived tidal tails. 

• (c) Apocentre: Gas density is increased in the cen¬ 
tre and suppressed in the outskirts. This produces a strong 
nuclear burst and mild off-nuclear star formation. 

• (d) Second Approach: Morphology is very similar to 
that at apocentre (c), but with lower levels of star formation. 

Only columns (b)-(d) correspond to the interacting 
phase (Section 12.311 . Merging and post-coalescence phases 
are omitted. The incoming phase (column a) is included 
for comparison (identical to the isolated case). It is evident 
that the spatial distribution of star formation in interacting 
galaxies has highly complex morphology. (Section [5] briefly 
discusses the ring-like feature in columns c and d.) 
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0 12 3 4 5 

time [Gyr] 


Figure 3. Global star formation rate (SFR) versus time (case 
study). Vertical lines delimit the stages of merging: first approach, 
interaction, coalescence, and post-coalescence. Red (blue) refers 
to the primary (secondary) galaxy, and black to the arithmetic 
sum of the two. Solid (dashed) curves indicate interacting (iso¬ 
lated) galaxies. Top: Global SFR. Middle: The logarithm of the 
SFR enhancement (SFR in interaction divided by SFR in iso¬ 
lation). Grey horizontal line indicates SFR enhancement equals 
unity. Bottom: Orbital separation. The galaxies experience two 
bursts, one in the interacting phase, the other at coalescence. 
The secondary galaxy exhibits weaker SFR and stronger SFR- 
enhancement (and suppression) than the primary. 

3.2 The Evolution of Global Star Formation 

Figure [3] shows SFR (top), SFR enhancement (middle), and 
separation (bottom) as a function of time (case study). SFR 
enhancement is defined as the SFR of the interacting 
galaxy divided by the SFR of its isolated counterpart. The 
solid vertical lines demarcate the different stages of merging 
(left-to-right): incoming, interacting, coalescing, and post¬ 
coalescence. The red (blue) curves refer to the primary (sec¬ 
ondary) galaxy, and the black curves represent their arith¬ 
metic sum (labelled ‘total’ in the Figure). Solid (dashed) 
lines represent interacting (isolated) galaxies. 

Left alone, the two galaxies in isolation experience a 
simple decaying star formation history. In interaction, on 
the other hand, the galaxies experience two bursts of star 
formation: one between first and second pericentric passage, 
and another when the final merger occurs. 

The Figure illustrates the nuisances described in Sec¬ 
tion m\ The spikes at second pericentric passage signal the 
overlapping of the two 10-kpc-radius spheres encompassing 



0 1 2 3 4 5 

time [Gyr] 


Figure 4. SFR versus time, in spherical shells of radii [0 — 0.3], 
[0.3 — 1], [1 — 3] and [1 — 10] kpc (light-to-dark blue), and in to¬ 
tal (black). Secondary galaxy, case study, interacting phase only. 
Solid (dashed) refers to the interacting (isolated) galaxy. Top: 
SFR per shell. Middle The logarithm of enhanced-SFR per shell. 
Grey horizontal line indicates SFR enhancement equals unity. 
Bottom: Orbital separation (dotted refers to snapshots not consid¬ 
ered). The innermost region dominates the total SFR, especially 
late in the interaction. Enhancement there is sustained, reaching 
factors higher than 10. The second shell is strongly enhanced for 
^1 Gyr, followed by suppression. The outermost shells are briefly 
enhanced, and quickly suppressed thereafter. 


each galaxy. For each galaxy, the SFR boost is caused by 
contamination from the ‘invading’ companion. This is par¬ 
ticularly dramatic during coalescence, as the two spheres 
merge, causing the red and blue curves to converge. Our 
aim here is to map star formation in the interacting phase, 
where the two distinct galaxies are clearly identified. For the 
rest of this paper, these nuisances are avoided and ignored. 

Our two metrics, SFR and SFR enhancement, are com¬ 
plementary. The larger galaxy has higher SFR and lower 
SFR-enhancement (compare red to blue). In other words, 
the smaller galaxy makes fewer stars, but is more suscep¬ 
tible to the encounter because of the relative tidal forces 
acting between the galaxies (e.g., iD’Onghia et al.ll201Cll ). No¬ 
tice that the interacting phase exhibits both enhancement 
and suppression. This is only true for configurations with 
sufiiciently-long interacting timescales. In general, this trend 
is interrupted by merging. 
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3.3 The Spatial Evolution of Star Formation 

Figure |4] shows the time evolution of SFR (top) and SFR 
enhancement (middle) in concentric spherical shells centred 
around the supermassive black hole of the secondary galaxy. 
We only discuss this galaxy for the sake of brevity. We focus 
on the following radii: [0 — 0.3], [0.3 — 1], [1 — 3], and [3 — 10] 
kpc (light-to-dark blue). The total SFR (within 10 kpc) is 
shown in black. Solid (dashed) curves refer to the interact¬ 
ing (isolated) case. The bottom panel shows the orbital sep¬ 
aration as a function of time. The dotted curve represents 
periods excluded from the analysis (outside the interacting 
phase and overlapping 10-kpc spheres, see Section E31) . 

The innermost region (with distance R < 0.3 kpc from 
the centre) exhibits the strongest levels of triggered star for¬ 
mation. It takes only ~0.5 Gyr after first pericentric passage 
for this region to account for half of the star formation in 
the galaxy. By ~1 Gyr after first passage, nearly all of the 
star formation is taking place in this region (compare light 
blue and black curves). This region is enhanced up to factors 
of ~15 (compare to the global SFR-enhancement of ~2-3, 
black curve). The second shell (0.3 < R < 1 kpc) experiences 
a slightly weaker burst - with shorter duration (~1 Gyr), 
followed by suppression. SFR in this region is enhanced by 
factors of >^5-6. The two outermost shells experience a brief 
(and weak) episode of SFR enhancement, and are quickly 
suppressed. In particular, the shell next to last (1 < R < 3 
kpc) experiences the strongest suppression - two orders of 
magnitude below the isolated case. 

In summary, the central regions experience stronger and 
longer periods of SFR enhancement, whilst activity in the 
outskirts is largely suppressed during the encounter. The in¬ 
crease in central star formation is due to the redistribution of 
gas is caused by non-axisymmetric tidal torques produced by 
the interaction, which leads to high gas density concentrated 
in the centre (Figure Hl), and which has been pre v iously seen 
in si mulations (e.g., iMihos fc Hernauisd 1 19961 : llono et ^ 
l2004h . This increase in nuc l ear SF R is i n line with recent sim - 
ulations bv iHopkins et all ll2013ll and iRenaud et ah] (l2015l l 
- although those analyses do not compare against SFR in 
isolated galaxies. To our knowledge, we are the first to re¬ 
port suppression of star formation at galacto-centric radii. 
In Moreno et al. (in prep) we investigate the mechanisms 
responsible for this suppression, and leave a discussion of 
that effect for this forth-coming paper. 


4 RESULTS: 75 MERGER SIMULATIONS 

4.1 Eccentricity & Impact Parameter 

We now explore various eccentricities and impact parame¬ 
ters: e = {0.85, 0.90, 0.95, 1.0, 1.05}, b = {2, 4, 8,12, 16} kpc. 


4.1.1 Global Star Formation Efficiency 

It is impractical to describe every possible merger in our 
suite with the same level of detail devoted to our case study 
(Section[3]). Instead, we adopt the integrated star forma¬ 
tion, which is equivalent to the total mass in stars created 
during the interacting phase (Section 12.311 : 

m..new=/ SFR{t)dt, (1) 

J interacting phase 


Secondary Goloxy (R < 10 kpc) 



0 1 2 3 4 5 6 7 8 91011121314151617 
Impoct Porometer [kpc] 


Figure 5. Global stellar mass, m*, new (equation^ top) and SFR 
efficiency, csfr (equation [2] bottom), for secondary galaxy in the 
interacting phase. Colours indicate eccentricities (purple-to-red): 
e = {0.85, 0.90, 0.95,1.0,1.05}. Vertical lines refer to these impact 
parameters: b = {2, 4, 8,12,16} kpc. Symbols are slightly offset for 
clarity. For all 25 orbits, more stars are made in the interacting 
galaxy than in its isolated counterpart (all symbols in the bottom 
panel are above the grey horizontal line). 


This q uantity is equivalent to the TSFR’ of iDi Matteo et ^ 
ll2007h . except that we integrate exclusively over the inter¬ 
acting period. Th is is also equivale nt to the mass in gas that 
turns into stars (IGox et al.l l2008ll . Alternatively, we could 
adopt SFRmax, the maximum SFR. However, the intermit¬ 
tent nature of star formation renders this quantity inade¬ 
quate because it depends strongly on our time-step choice. 

We also introduce the star formation efficiency: 


esFR = 


( 2 ) 


where m*_new is defined in equation 0, and is its 

analogue in i solation. T h is qua ntity is similar to the ‘burst 
efficiency’ of IGox et al.l ll2008ll - except that (1) we only 
consider the interacting phase; (2) they compute a difference 
where we compute a ratio; and (3) they use the total SFR of 
the two galaxies, whilst we use it for individual galaxies (this 
section), and for subregions therein (Section 14.1.21 below!. 

Figure [5] shows the integrated star formation (top) and 
star formation efficiency (bottom). The solid vertical lines 
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Secondary Galaxy (R < 1 kpc) 


Secondary Galaxy (1 < R < 10 kpc) 
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Figure 6. Stellar mass {top), SFR fraction (equation (3] middle) and SFR efficiency (bottom) in the secondary galaxy: nuclear (R < 1 
kpc) versus off-nuclear (1 < R < 10 kpc) contributions (left versus right). Colours, symbols, and vertical lines as in Figure [5] Grey 
horizontal line refers to SFR efficiency of unity. Star formation is split into ~87 % (nucleus) and ~13 % (outskirts). SFR efficiency is 
enhanced (suppressed) inside (outside) the central kpc. 


mark impact parameters (the symbols are offset for clarity), 
and the colours (indicated in the key) represent eccentrici¬ 
ties. All of our 25 orbits produce m*,new ~ (2 —6) x 10® M©, 
with efficiencies ranging from ~1.5—3. 

It is beyond the scope of this work to identify exactly 
how star formation tr iggering depends on e and b (see, e.g., 
iDi Matteo et al.]l2007l ). In broad terms, the amount of star 
formation is governed by two factors: the strength of the 
interaction and its duration. In particular, highly eccentric 
orbits with large impact parameters last longer, leading to 
larger values of m*,new (top panel). Correcting for the du¬ 
ration of the orbit (by dividing by shows that low 

eccentricities and intermediate impact parameters lead to 
the highest star-formation efficiencies (bottom panel). 

4.b.3 Nuclear versus Off-Nuclear Efficiency 

We split our galaxies into two parts: the nuclear region 
(R < 1 kpc) and the off-nuclear region (1 < R < 10 kpc). 
Figure [S] is analogous to Figure [S] but constrained to the 
nuclear (left) and off-nuclear (right panels) regions. 

The fraction of stellar mass per region, defined as 

_ /region P*’ new(R) dR 

/o°‘'‘’"47rR2p*.new(R)dR’ 


is shown in the middle panels. p*,new(R) is the mass-density 
radial profile of new stars created throughout the interac¬ 
tion. Our 25 mergers produce m*,new ~ (1.5 - 5) X 10®M© 
in the central kpc (upper left), and ~ (0.5 — 1) x 10®M© in 
the outskirts (upper right). Across all orbits (all values of e 
and b), the mass in new stars is ~87% in the nucleus (mid¬ 
dle left), and ~13% elsewhere (middle right). These portions 
are nearly independent of e and b (middle panels). 

In stark contrast to the global case (Figure [5|), star 
formation efficiency is not above unity across the entire 
galaxy. Instead, it ranges between ~4-10 in the nucleus and 
~0.2—0.5 in the outskirts. In other words, interactions en¬ 
hance star formation in the centre, and suppress it at large 
radii. In particular, the suppression found in our case study 
(Figure |4|| is generic across all of our 25 merger simulations. 


4.2 Comparing the Two Galaxies 

Figures [5] and [6] demonstrate that quantities like m*,new, 
fsFR, and esFR are weakly dependent on e and b. Hereafter 
we only report averages over our 25 simulations (denoted 
in brackets). Error bars represent standard deviation of the 
mean. We now include the primary galaxy in our analysis. 
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Figure 7. Mean mass in new stars, (m*,new) {iop)-) mean SFR 
fraction (fspR)? equation [5] middle), and mean SFR efficiency, 
(csfr) (equation [JJ bottom)- across 25 merger simulations (in¬ 
teracting phase only), with varying e (eccentricity) and b (im¬ 
pact parameter). We only include orbits in the “e” orientation. 
Blue (red) lines correspond to the secondary (primary) galaxy. 
Grey horizontal line refers to mean SFR efficiency of unity. Solid 
vertical lines indicate the following regions (concentric spherical 
shells): [0 — 0.3], [0.3 — 1], [1 — 3] and [3 — 10] kpc. Symbols are dis¬ 
placed for clarity. Both galaxies have centrally-concentrated star 
formation, enhanced efficiencies in the central regions (first two 
bins), and suppressed efficiencies (below unity) everywhere else. 
Trends are more evident for the secondary galaxy. 



Figure 8. Mean star formation efficiency, (csfr) (equation [4|l. 
across 75 merger simulations (interaction stages only). Blue (red) 
lines represent the secondary (primary) galaxy. Grey horizontal 
line refers to mean SFR efficiency of unity. Vertical lines indi¬ 
cate the following regions (concentric spherical shells): [0 — 0.3], 
[0.3—1], [1 — 3] and [3 — 10] kpc. Symbols are displaced for clarity. 
Three orientations are considered (25 orbits each): “e” (aligned), 
“f’ (perpendicular), “k” (anti-aligned). See Figure ^ and Ta¬ 
ble ^ for definitions. Both the “e” and “f’ produce centrally- 
concentrated star formation in the secondary galaxy. Only the 
“e” orientation has this effect on the primary. For anti-aligned 
galaxies (“k” orientation), the effect on either galaxy is minimal. 

mark the four spatial regions of interest, and values are off¬ 
set for clarity. For both galaxies, the presence of a compan¬ 
ion makes star formation more centrally concentrated. The 
bottom panel shows that (esFu) is enhanced in the inner re¬ 
gions, and suppressed in the outer shells. This is particularly 
evident for the secondary galaxy (by a factor of ~0.3 in the 
last shell), and less obvious for the primary galaxy (by a fac¬ 
tor of ^0.7 in that same shell), suggesting that the smaller 
galaxy is more susceptible to interaction-driven effects. 


We define the mean star-formation efficiency as 

M (4) 

\''t'*, new/ 

and the mean star-formation fraction as 


^/region P*. new(f?) dR) 

P*,new(i?) dR)' 


(5) 


(Strictly speaking, (fspn) and (esFa) are not proper aver¬ 
ages, but ratios of averaged quantities.) We concentrate on 
four spherical shells spanning the following radii: [0 — 0.3], 
[0.3 — 1], [1 — 3] and [3 — 10] kpc (as in Figure [4]). 

Figure [7] shows (mt,new) (upper panel), (Isfr) (middle 
panel) and (bsfr) (bottom panel). Blue (red) lines repre¬ 
sent the secondary (primary) galaxy - solid (dashed) lines 
correspond to interacting (isolated) galaxies. Vertical lines 


4.3 Disc Spin Orientation 

Figure [S] shows the mean star formation efficiency (equa- 
tion|4)| for the “e” (solid lines, aligned discs), “f’ (dot-dashed 
lines) and “k” (dotted lines) orientations. These choices rep¬ 
resent (nearly) aligned, perpendicular and anti-aligned discs 
(Figure [1] and Table [T]). The solid lines here are identical to 
the bottom panel of Figure [7] 

For the secondary galaxy (blue), the “e” and “f’ ori¬ 
entations behave similarly: star formation efficiency is en¬ 
hanced in the centre, and suppressed in the outskirts. In 
contrast, the “k” orientation exhibits substantially weaker 
enhancement (suppression) in the nucleus (outskirts). For 
the primary galaxy (red), the “e” orientation is the only 
one showing strong enhanced (suppressed) star formation 
efficiency in the nucleus (outskirts). The other two orienta¬ 
tions, are weakly enhanced in the central regions. The “f’ 
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orientation is also weakly enhanced in the outskirts, whilst 
the “k” orientation is consistent with unity in that region. 

The “k” orientation (both galaxies) and “f’ orientation 
(primary galaxy) exhibit the weakest deviations from the 
isolated setting. After checking individual cases, these devi¬ 
ations are caused primarily by orbits with h = 1,2, which 
allow disc interprenetration at first pericentric passage. 


5 DISCUSSION & FUTURE WORK 
5.1 Emerging Picture 

Standard wisdom suggests that g alaxy enc o unters 
trigger bursts of star forma tion (iHernauistI 1 19891 : 
iBarnes fc HernauistI Il99ll . Il996li . However, the spatial 
extent of such episodes, and whether or not these events oc¬ 
cur in both galaxies, remain poorly understood. To address 
these issues, we employ a suite of 75 merger simulations. 
We encapsulate the level of star formation in each orbital 
configuration in terms of the star-forming efficiency, defined 
as the ratio of the mass created during the interaction 
and the equivalent amount created in isolation ((espn), 
equation Hi. We find that, for the entire galaxy (Figure O 
and for individual shells (Figure |3, this quantity is weakly 
dependent on t and 6 - at least for the range of values 
explored he re (selected to be consiste nt with cosmological 
simulations, iKhochfar fc Burkertll20()^ . 

The intensity and spatial extent of star formation in 
interacting galaxies is strongly driven by the alignment 
of the two disc spins (Figure [H]). With the exception of 
a few anomalous cases where the two galaxies interpene¬ 
trate one another, interactions with anti-aligned spins (“k”- 
orientation) have a minimal effect on the participating 
galaxies. The alignment of the two spins, however, per¬ 
mits centrally-concentrated star formation. The secondary 
galaxy experiences nuclear starbursts with either perpen¬ 
dicular (“f’-orientation) or aligned (“e”-orientation) spins. 
The primary galaxy, on the other hand, requires the two 
spins to be aligned for nuclear star formation to occur. 

Our results suggest that relative spin disc orientation is 
the dominant factor behind centrally-concentrated star for¬ 
mation in interacting galaxies, particularly for the smaller 
galaxy. Interestingly, we find that whenever enhanced star 
formation in the nucleus is triggered, this is always accompa¬ 
nied by suppression of activity at large galacto-centric radii. 
The detailed physics behind this process is explored in forth¬ 
coming work (Moreno et ah, in prep). 


5.2 Connection to Other Work 


Little systematic work has been done to quantify the spatial 
distribution of star formation in mer gers (i.e., by exploring a 
broad range of orbital parameters). Di Matteo et al.l (l2007l f 
is an exception. Unfortunately, these authors only report 
maps at a few time snapshots, with limited quantitative in¬ 
formation on how intense this process is at each region. 

Interestingly, a large number of observational and the¬ 
oretical works emphasize the importance of widespread star 
formation - both in tidal tails and bridges. These include 


the Antennae Ga l axies NGC 4038/39 (lMiTabel_et_^ 


I Wang et al.l l2004l : iRenaud et al.l l2008l : Brandi et al.l 


1998 


200C 


Karl et al. I I 2 OIOI: iTevssier et al. Il2010l: IWhitmore et al.ll2010l: 


Zhan^^^^EwJjKa^^^L gOU . 2^31: Privon et ahlboiS 


Renaud_et_al .ll2014l:lRenaud et al.feoi^: the Mice 4676A/B 

dSotnikova fc Reshetnikov _ W9^ Read I 2 OO 3 I : iBarned l2004l : 


IPrivori et al.l 20]T Wjld_e^^ 201^: An 299 1C 69 4./NG C 


3690 I Alonso-Herrero et al. 2000l. 20091: ISl iwa et al.l |201^: 


NGC 151 2/151 0 ( Koribalski fc Lopez- Sanche^ 200^ : 
Ducd et al.l 2014 |l: N G C 2201 / IC 2163 (j Elmegre en et al.l 


2 OOOI: Struck et al.| l2005l: lElmegreen et ^ I 2 OO 6 I : 

Kaufman et al.l 2012_|Mineo_et^L 2013ll : and the Con¬ 
dor NGC 6812 ( Eufrasio et al.l 20l4) . Whilst compelling. 


it is not clear if the lessons learnt from these specific cases 
can be generalised. 

One could argue that tail/bridge star formation is ubiq¬ 
uitous. In our simulations, this only occurs for orbits where 
the two discs interpenetrate, and only for a few Myr. Our 
analysis suggests that this mode of SFR localization is sub- 
dominant (Figure [H| l. This is consistent with observations. 
ISchmidt et al.l (l2013h embarked on a systematic study of 60 
merging systems (3d-hst). They show that only ~3% of 
their mergers exhibit star formation in the regions “in be¬ 
tween” - showing that, whilst this phenomenon does exist, 
it is rather uncommon. For this reason, we do not attempt 
to refine our underlying modelling, nor do we try to compare 
against any of the above sp ecialised studie s. _ 

It is worth noting that ISchmidt et al.l ll2013ll also com- 
pare to simulations, wh ich happen to be very similar to ours 
JCox et al.ll2006l . I 2 OO 8 II . However, they claim that their sim¬ 
ulations do not agree with observations. They report more 
instances with both galaxies exhibiting elevated levels of star 
formation simultaneously in the simulations than in their 
observations (~59% versus ~32%). They argue that this is 
due to the fact that both galaxies in their simulations have 
identical initial gas fractions; which might not be the case 
for real interacting galaxies. This interpretation is certainly 
viable. Another explanation is that their simulations and ob¬ 
servations could be covering a non-representative handful of 
time snapshots and projections (observational viewpoints) - 
which might be incapable of fully capturing the interaction- 
induced star-forming process in these systems (60 observed 
and 296 simulated maps). In our simulations we see situa¬ 
tions where both galaxies begin with the same gas fraction, 
but exhibit different levels of star formation at a given point 
in time. For instance, in our case study, the primary galaxy 
is still forming stars during the late stages of the interacting 
phase, whilst the companion is not (upper panel of Figure|3]). 


5.3 Future Directions 


There are several directions worth pursuing. Figure [2] shows 
the complex morphology of star forming gas in interacting 
galaxies. In particular, the ring-like structure appearing af¬ 
ter first passage is intriguing. If robust, this imprint could 
provide a promising method for identifying galaxies that 
have experienced a close encounter in the past - e.g., com¬ 
pare columns (c) and (d) in Figure [2] to column (a). This 
could prove to be a powerful alternative to selecting interac¬ 


tions via the presence of tidal tail s (|KiMdaItepe_et_aJj 


Dare et ^ l2010al : ICasteels et al.l I 2 OI 3 I : Hung et ^ 


2012 


20131. 


2014f) . which dissipate faster than this ring-like feature (col¬ 
umn b) - provided that observations are capable of measur¬ 
ing these detailed features (via Ha and [OHI] emission-line 
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maps). We acknowledge that at this point, it is not clear 
if the properties of this feature (or its existence and ob¬ 
servability) depend strongly on numerical/modeling effects. 
This warrants a more rigorous study. The conclusions of this 
paper are not affected by this ring, which only accounts for 
~4% of the produced stars in our simulations. 

We also plan to co n struct a ‘mock survey’ - along the 
lines of ISchmidt et al.l ll2013h . but with more snapshots, 
viewpoints, and merger configurations. This framework will 
replace three-dimensional spherical shells with cylindrical 
annuli, to compare better with observations. Azimuthal de¬ 
pendences within these annuli will capture tidal tails and 
other features. This mock survey will be better suited for 
capturing the localization of SFR on a snapshot-by-snapshot 
basis - thereby overcoming some of the limitations inherent 
to the time-integrated metrics (e.g., csfr, equation [2} em¬ 
ployed in this paper. 

Our results allude to the importance of exploring other 
orientations and mass ratios. We find that the smaller 
galaxy has a more evident (nuclear enhanced, off-nuclear 
suppressed) response to the encounter. With this, we spec¬ 
ulate that this trend will continue for other mass ratios (be¬ 
yond our adopted choice, ~2.5:1). That is, for more dis¬ 
crepant ratios, nuclear triggering in smaller galaxy might 
require weaker alignments, and stronger alignments for the 
larger galaxy. Clearly, this requires exploring a more refined 
range of orientations, and extending our mass-ratio regime. 

Ideally, it is also desirable to explore other mass, gas 
fraction, and bulge-to-disc ratio combinations. We must em¬ 
ploy extra care becaus e massive galaxies tend to be more 
bulge -dominated (e.g., iKauffmann et aO 2QQ^Ku^^t^^] 


2014|) , and have lower gas fractions (e.g.. lStewart et ^120091 : 


Catinella et al. I [2^. Furthermore, bulge- dominated sys¬ 


tems tend to inhabit de nser environments l|Dressleil Il980l : 


merg 

er” approximation less valid 

Martie & BournaudI 

boos 

: lMorencll2012l: iMoreno et al.ll2013 

)• 


6 CONCLUSIONS 

The aim of this work is to investigate the spatial extent 
of star formation in interacting galaxies. To address this, 
we employ a suite of 75 idealised SPH merger simulations 
(with stellar masses M* = 14 x 10^°Mq and 5.7 x 10® Mq, 
bulge-to-disc mass equal to 0.31, and gas fraction equal to 
0.25). We only consider three relative spin orientations, la¬ 
belled “e”, “f ’ and “k”. These choices are meant to represent 
two nearly aligned discs (“e”-orientation), two perpendicu¬ 
lar discs (“f’-orientation), and two anti-aligned discs (“k”- 
orientation). For each fixed spin disc orientation, we explore 
five eccentricities (e = {0.85,0.90,0.95,1.0,1.05}) and five 
impact parameters (6 = {2, 4, 8, 12,16} kpc). 

In this work, we quantify the spatial distribution of 
star formation in terms of spherical shells of radii [0 — 0.3], 
[0.3 — 1], [1 — 3] and ]3 — 10] kpc. We focus exclusively on 
the interacting stage, between first and second pericentric 
passage. We employ star-formation efficiency (equation [2]) 
- defined as the integrated star formation during the 
interacting phase, divided by its analogue in isolation - to 
encapsulate interaction-induced effects for each orbit. We 
warn that our results only apply to the period when the two 


merging galaxies can still be identified as separate entities. 

Our main results are: 

• Interactions generally produced enhanced star forma¬ 
tion in the centre, and suppressed activity in the outskirts. 
This effect is weakly dependent on the values of e and b 
probed here (Figures [SJ. 

• These trends are more pronounced in interactions with 
strongly aligned disc spin orientations (Figure |8}, particu¬ 
larly for the secondary (smaller) galaxy (Figure [T]). 

It is our hope that the numerical investigations pre¬ 
sented in this paper motivate detailed observational stud¬ 
ies with samples of galaxies drawn from ongoing and future 
deep-held and integral-held-spectroscopic surveys. 
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